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We have previously shown that quantum-enhanced atom interferometry can be achieved by map¬ 
ping the quantum state of squeezed optical vacuum to one of the atomic inputs via a beamsplitter-like 
process [Phys. Rev. A 90 , 063630 (2014)]. Here we ask the question: is a better phase sensitivity 
possible if the quantum state transfer (QST) is described by a three-mode-mixing model, rather 
than a beamsplitter? The answer is yes, but only if the portion of the optical state not transferred 
to the atoms is incorporated via information recycling. Surprisingly, our scheme gives a better 
sensitivity for lower QST efficiencies, and with a sufficiently large degree of squeezing can attain 
near-Heisenberg-limited sensitivities for arbitrarily small QST efficiencies. Furthermore, we use the 
quantum Fisher information to demonstrate the near-optimality of our scheme. 
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I. INTRODUCTION 

Squeezed light is a resource that enables the detection 
of phase shifts below the standard quantum limit (SQL) 
1/V^, where N is the number of detected photons [1]. 
Famously, it was recently used to enhance the sensitiv¬ 
ity of the GEO 600 gravitational-wave detector [5] and 
the Hanford LIGO detector [2; this speaks to the ma¬ 
turity of squeezed-light-generation technology. As a con¬ 
sequence, squeezed light is a promising controllable re¬ 
source for generating squeezed atomic states - and there¬ 
fore enabling sub-SQL atom interferometry - via a quan¬ 
tum state transfer (QST) process that maps the state of 
the squeezed light to the atomic field mil]. Given the 
low fluxes of current atomic sources (compared with pho¬ 
ton sources), and the technical barriers to increasing this 
flux [T^ IT5] , it is likely that atom interferometers could 
only measure potential violations of the weak equivalence 
principle |14H16j . quantum gravitational effects HZ], and 
gravitational waves da [IS] if practical operation below 
the SQL is achieved. 

Unfortunately, in practice QST between atoms and 
light is imperfect; any transmitted component of the op¬ 
tical squeezed state behaves as a loss mechanism that 
can drastically reduce the effectiveness of the squeezing 
da HI]. However, the deleterious effects of imperfect 
QST can be reduced if the transmitted photons are in¬ 
stead measured, and the information obtained is com¬ 
bined in just the right way with the atom-interferometer 
signal. This technique is called information recycling 
[22j . and we have previously shown its effectiveness in 
atom interferometry schemes enhanced by either single¬ 
mode or two-mode squeezed optical vacuum [53] , and also 
for a general class of Heisenberg-limited input states [23] . 
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The majority of this previous work assumed that the 
QST process behaved as an atom-light beamsplitter with 
a beamsplitter reflection coefficient Q corresponding to 
the efficiency of the QST process (e.g., when Q = 1 all 
photons are mapped to atoms, when Q = 1/2 only 50% of 
the photons are mapped to atoms, etc.). This description 
is only valid when the mean number of photons in the in¬ 
put optical state is much smaller than the total number 
of atoms. In contrast, when the number of photons and 
atoms are comparable, the QST process is more appro¬ 
priately described as a three-mode mixing process (with 
one photon mode and two atomic modes corresponding 
to the mode initially occupied and the mode into which 
atoms are outcoupled during QST). We have previously 
investigated this QST process in the regime of moderate 
squeezing, where there is a quantitative but not qualita¬ 
tive difference to a beamsplitter QST process [23] . 

In this paper, we consider this three-mode mixing 
regime in considerably more detail, including the large 
squeezing regimes where the number of squeezed pho¬ 
tons is comparable to the number of atoms, yielding dy¬ 
namics and sensitivities that are both quantitatively and 
qualitatively different to the beamsplitter model of QST. 
In particular, we show that if a single-mode squeezed 
optical vacuum is partially mapped to two atomic in¬ 
puts of a Mach-Zehnder (MZ) interferometer, and the 
transmitted photons are monitored and information re¬ 
cycling is applied, then lower QST efficiencies can in fact 
lead to better phase sensitivities. Indeed, for a suffi¬ 
ciently large squeezing parameter, we show that near- 
Heisenberg-limited phase sensitivities are obtainable for 
arbitrarily low QST efficiencies. Crucially, this result is 
somewhat different to [25], which showed that a MZ in¬ 
terferometer can achieve Heisenberg-limited sensitivities 
if the two inputs are a coherent state and a single-mode 
squeezed vacuum state. In our scheme, the three-mode 
mixing QST process results in some non-trivial interfer¬ 
ometer input state that is entangled with the transmitted 
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photons. We also demonstrate the near-optimality of our 
phase-estimation procedure by computing the quantum 
Fisher information. 


II. THEORETICAL MODEL 
A. Atom-light QST process 

The simplified atom-light coupling model used 
throughout this paper is derived in considerable detail 
in [53] (see also EHin]), so we only briefly summarize 
it here. We consider an interactionless Bose-Einstein 
condensate (BEC) [26] comprised of atoms with two 
hyperfine ground states |1) and |2). These two states 
are coupled with a Raman transition via two counter- 
propagating optical fields, detuned from some excited 
state |e) that is not populated on the timescale of this 
coupling (see Fig. [^. The optical field addressing the 
12) o |e) transition is assumed to be a bright coherent 
state, and is therefore treated classically. Energy and 
momentum conservation ensures that the Raman tran¬ 
sition has a high degree of mode selectivity. So, given 
a BEC that initially occupies a single motional mode of 
state |1) (mode di), during a Raman transition predom¬ 
inantly one mode of the electromagnetic field addressing 
the |1) —>■ |2) transition (mode b) outcouples atoms to a 
different motional state in |2) (mode 02 ). Here di and 
b are bosonic modes obeying the usual commutation re¬ 
lations. The effective interaction Hamiltonian describing 
this process is 


H = hg ^did^d -I- H.c.^ , (1) 

where g is the coupling strength between light and atoms, 
and we have assumed that the system is on-resonance 
such that the energy transferred by the two-photon tran¬ 
sition perfectly matches the change in electronic and ki¬ 
netic energies of the atom. 

Typically, the occupation of modes d 2 and b is assumed 
to be much less than the occupation of mode di. In this 
case we can treat di as an undepletable reservoir and 
make the approximation di —>■ where Na^ is the 

mean number of atoms in mode di. Consequently, Eq. Q 
reduces to the beamsplitter Hamiltonian 

Hbs = (dj& + i?.c.) , (2) 

which corresponds to the linear evolution 

d 2 (t) = - Q{t) a2{to) - i\JQ{t) b{to), (3a) 

d(t) = - *V^Q(i)d2(to) • (3b) 

Here Q{t) = sin^{g is interpreted as the reflection 
coefficient of this atom-light beamsplitter, and therefore 
quantifies the efficiency of QST between the photon mode 
b and atomic mode di. For example, if Q{t) = 1 then 



FIG. 1. (Color online) Energy level scheme for a three-level 
Raman transition comprising two non-degenerate hyperfine 
ground states, |1) and |2). The BEC is initially formed in the 
state |1), and population is transferred to |2) via the absorp¬ 
tion of a photon from mode b and the emission of a photon into 
the optical field addressing the |e) |2) transition (which we 

treat as a large classical field). Since the detuning A of the 
optical fields is large, the excited state |e) is not appreciably 
populated on the timescale of the QST dynamics. 


the quantum state of b(to) is completely mapped to d(t), 
whereas if Q{t) = 0.5 only 50% of the photons couple to 
the atoms via the Raman transition. 

When the number of photons is comparable with the 
number of atoms in the condensate, the undepleted reser¬ 
voir approximation cannot be applied, and the QST pro¬ 
cess must be described by the three-mode mixing dynam¬ 
ics 

iai = ga 2 b\ (4a) 

ia2 = gaib, (4b) 

ib = ga 2 a\, (4c) 


which are the Heisenberg equations of motion assuming 
Hamiltonian Q. Although an analytic solutio n is po ssi- 
ble for times much shorter than g~^ (see Sec. IVA), in 
general Eqs. @ must be solved numerically. 

Although beamsplitter dynamics no longer apply, the 
QST process can still be quantified via the generalized 
QST efficiency 


^ {d\{f)d2{t)) 

(6t(tp)6(to))’ 


(5) 


which compares the number of atoms outcoupled com¬ 
pared with the total number of input photons. 

This three-mode mixing process Q has been pre¬ 
viously studied within the context of non-degenerate 
parametric down-conversion [23 [28] , and recently by us 
within the context of squeezed-light-enhanced atom in¬ 
terferometry [23] . This prior analysis of three-mode mix¬ 
ing QST dynamics was in regimes which are qualitatively 
similar to a beamsplitter QST process. Below we consid¬ 
erably extend the analysis of the single-mode squeezed- 
light-enhanced atom interferometer presented in [23] : in 
particular, we investigate the regime of moderate to large 
optical squeezing, where the three-mode mixing QST dy¬ 
namics is qualitatively different to a beamsplitting QST 
process. 
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FIG. 2. (Color online) A hybrid atomic-photonic interferometer enhanced with single-mode squeezed optical vacnum and 
information recycling. Left panel: Atoms in mode Si are (partially) outconpled to mode 02 via a stimu lated Raman transition, 
with the goal of mapping the quantnm state of the squeezed optical vacuum b{to) to a 2 {ti) (see Sec. IIAl. Right panel: the 
two atomic modes Si and S 2 form the two arms of a MZ atom interferometer; any relative phase shift (j) that accrues between 
the two modes during the evolution of the atoms through the interferometer can be determined by measuring the number 
difference, Sa, of the two atomic modes at time tf. By interfering any light transmitted after the QST process [fe(ti)] with a 
local oscillator (L.O.), assumed to be a large coherent state, a measurement of the phase quadrature Yf, of b{ti) is achieved. 
Information recycling is then incorporated as a post-processing step, whereby the information-recycled signal S is constructed 
from the atomic and photonic signals, Sa and Sb respectively. The signal 5 generally gives a more sensitive phase measurement 
than the purely atomic signal Sa- 


B. Hybrid atomic-photonic interferometer with 
information recycling 


The atom-light QST process and information recycling 
is incorporated into a MZ interferometer as shown in 
Fig. i Initially (at time to), all the condensate atoms 
are in a coherent state in mode di, and our photon mode 
is assumed to be a single-mode squeezed vacuum state 

m 


d(to) = i?coshr -I- sinhr, 


( 6 ) 


where d is a vacuum mode and r > 0 is the squeezing pa¬ 
rameter. Our state is squeezed in the phase quadrature, 
as this is optimal for our scheme. Using an appropri¬ 
ately chosen coherent optical held, the optical mode b is 
coupled to modes di and 0,2 via a Raman transition, as 
described in Sec. |II A with the aim of mapping the state 
of b{to) to 02 (^ 1 ) (see left panel of Fig.[^. 

After the QST process is complete, the two atomic 
modes di(ti) and 02 ( 11 ) form the inputs of a MZ atom 
interferometer (see right panel of Fig. [^. Here light 
pulses are used to coherently interfere, redirect, and re¬ 
combine the two atomic modes. The evolution of the 
atomic modes through the interferometer is most con¬ 
veniently described in terms of the pseudo-spin angular 


momentum operators: 

Jx = lialai + d|a2), 

Jy = - a\a2), 

Jz = - 0202). 


(7a) 

(7b) 

(7c) 


In particular, if the two modes accrued a relative phase 
shift (j) during periods of free evolution within the inter¬ 
ferometer, then this phase shift can be determined by 
measuring the relative number difference at time tf: 


Sa = 2Jxitf)=2 


J,,{ti)cos(j) 


Jx(ti) sint/) 


( 8 ) 


If there is complete QST, Q = 1, then the minimum 
sensitivity of this phase measurement can surpass the 
SQL [I1I13]: 


A(j) = 


| 50 ( 5 ,)| 




(9) 


where Nt is the total number of atoms in the initial BEG, 
V(O) = (O^) — (O)^ is the variance, and we have assumed 
that Nh = {b\to)b{to)) = sinh^ r ^ Nt- 

However, for the typical case of Q < 1, the squeezed 
optical state is only partially mapped to 0 - 2 - Since, to 
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first order, the QST process is a beamsplitter that mixes 
vacuum [ 02 (to)] with the squeezed optical input &(to) (i.e. 
a linear loss mechanism), the effects of the squeezing are 
drastically reduced, and the phase sensitivity is degraded. 
Fortunately, the lost portion of the metrologically use¬ 
ful quantum correlations still exist in the mode of the 
transmitted light, b{ti). As discussed in (53], if the phase 
quadrature of the transmitted light Yb = *[^(ii)—&^(ii)] is 
measured via homodyne detection, then the information- 
recycled signal 

s = VQSa - VT^Sb, ( 10 ) 

where Sb = \l{Na^ {ti))Yb, significantly ameliorates the 
effect of incomplete QST on the phase sensitivity. Again, 
the benefits of information recycling on this hybrid atom- 
photonic interferometer were only quantified for a) a 
beamsplitter QST process in the regime sinh^ r A^t 
and b) a three-mode mixing QST process in the regime 
of moderately-sized r where the QST process behaves 
qualitatively as a beamsplitter. Below, we consider in 
detail the regime where Nb > Nt, and show that infor¬ 
mation recycling gives surprisingly good sub-SQL phase 
sensitivities with a counterintuitive dependence on r and 

Q. 



Q 


FIG. 3. (Color online) The QST-dependence of phase sen¬ 
sitivity Eq. ( |12[ ), which assumes a beamsplitter QST pro¬ 
cess, for different squeezing factors r and an initial con¬ 
densate of Nt = 10^ atoms. The minimum phase sensi¬ 
tivity of is only attained if r > Ccrit and 

Q « 2exp(-2r)-v/iVt- 


III. BEAMSPLITTER QST PROCESS 


Before considering the full three-mode mixing QST dy¬ 
namics, we first consider the simplest case where the QST 
process is a beamsplitter [see Eqs. ([^]. Here, the phase 
sensitivity of the information-recycled signal S can be 
computed analytically; although the beamsplitter QST 
model is invalid except in the regime where sinh^ r Nt, 
this provides some basic intuition and develops a base¬ 
line to compare with the more counterintuitive results 
of Sec. |IV[ Specifically, at the optimal operating point 
(j) = 7r/2, the variance and slope of the signal S are 


E(5) = Nte-"^^ + Q{Q - ginh^ r, (11a) 
d 4 ,{S) = VQ{Nt - 2Qsinh^r), (11b) 

respectively, yielding a phase sensitivity of [23] 


A())bs(g Q) 


+ {Q- e-^^) sinh^ r 
Nt — 2Q sinh^ r 


( 12 ) 


In the limit sinh^ r <C Nt, ~ exp(—r)/-\/Q7Vt, 

which surpasses the SQL provided Q > exp(—2r). 
However, in the intermediate regime where r is large 
enough that Q exp(—27-), and the second term un¬ 
der the square root in Eq. must be kept, but Nt » 
2Q sinh^ r: 

\/%^ + 2sinhV 

A<(.Bs(r, Q) « ^ -• (13) 


Equation (13) attains a minimum 


of 


A</>Sr(2) 


Nt-^VNt + ^ 


N, 


5/4 


N, 


3/4 ’ 


(14) 


at an optimal squeezing parameter of r°P‘(Q) « 

ln(4iVt/Q2)/4. 

We can therefore identify two qualitatively different 
regimes, delineated by rcrit = ’'bs*( 2 = !)■ When 
r < Tcrit, the minimum sensitivity only occurs at Q = 1, 
and is bounded by < A</)gg" < exp(—r)/i/]Vt; 

generally this upper bound is a good approximation to 
A<(’bs'- When r > Ccrit a minimum phase sensitivity of 
A^gg" « is always possible, however this occurs at 

a QST fraction less than unity (see Fig.[^. That the min¬ 
imum phase sensitivity requires a reduction in QST effi¬ 
ciency, and that furthermore Qopt ~ 2 exp(—2r)-\/]Vt —>• 0 
as r gets large, is a curious feature of our hybrid atomic- 
photonic interferometer. However, for a beamsplitter 
QST process an imperfect QST efficiency does not give 
an improved sensitivity over the sensitivity at Q = 1 and 
r = r^rit. Consequently, there is no inherent advantage in 
operating in a small Q, large r regime for a beamsplitter 
QST process. This is in contrast to a three-mode mixing 
QST process, which, as we show below, can achieve a 
near-optimal Heisenberg-limited phase sensitivity in the 
small-Q regime. 
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IV. THREE-MODE MIXING QST PROCESS 
A. Analytical solution in small-QST regime 

Although the dynamics of the three-mode mixing QST 
process cannot be analytically solved in general, an ap¬ 
proximate solution does exist for times ti ^ g~^ (equiv¬ 
alently, for Q ^ 1). Defining t = gti, we can apply the 
improved Euler (i.e. Heun’s) method [3D] to Eqs. @ and 
obtain solutions valid to second order in r: 

ai(r) = 0 , 1 - iTpa 2 + ^ai{Na^ - N^) + 0 {t^), (15a) 

a2(r) =02- irboi - ^a2(fVoi+ iVf, -|- 1) -|- 0{t^), (15b) 
S(r) = b-iTd\a 2 + ^b{Na 2 - -^aj + (15c) 


For notational compactness, we have written ai = ai(0), 
etc., Nai = d\di, and = b^b. Note that an expansion 
to second order in r is required in order to accommodate 
a non-zero QST efficiency, since Q = Ntr"^ + 0 {t^). 

Noting that the average number of particles in mode 
hi at time r is 

A^ai(r) = iVt(l-T^[l-^ (3cosh(2r) -|- l)]sinh^r), (16) 

at the optimal operating point (j) = 7r/2, the derivative 
and variance of the information-recycled signal S are 

d,f,{S) = \fWtT {Nai (r) - sinh^ r) (17) 

and 


ViS) = e-2"(l - NtT^)Na,iT) [l + NtT^ - 1) - l)] - NtT^e-’'sinhr^Nt{l - NtT^)NaAT) 

X {4 - ^ (4At -h 3e-^^ - 1) + fs - 2) -f coshr -t 9e-'^^ - 7] } 

-I- — 2NtT^e~^ sinhr + x -b At (3 -f 2 (l — 6e“^’') sinh^ ^) + § sinh^(2r)] 

sinh^ r [2 -b A^ -h (3e-’’ (3 - -b e’’ (36^’' -b ll)) -b 6cosh(2r)] 

sinh^ r [2Aj^ (3cosh(2r) -b 1) -b At (15 [2 cosh(2r) -b cosh(4r)] -b 11) -b | cosh^ r (35cosh(4r) -b 13)] |, (18) 
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respectively. Substituting these quantities into the linear error propagation formula Eq. and performing a power 
series expansion around r = 0, gives the following expression for the phase sensitivity, valid at small times r: 

— p(q 4 p ll 

A( 5 iQ^i(r, t) ="I 1 -b -T^ sinh^ r-b ^ At (At-b 3 — 2e“^’’)-b 2 (5 cosh(2r) — 6) sinh^ r ^+ 0 { t ^). (19) 


NtT\ 


In the r > 1 regime, the first, second and third terms 
in Eq. (19) scale as e~'~ jr, re’’ and {re^Y, respectively. 
A(/) will therefore approach infinity as r —>■ 0 unless we 
choose optimal r = = ln(C/r) for some constant C, 

in which case we obtain a minimum sensitivity of 




X/^3 I "ip I l-tQ^/8 

32^ -I- 8^ ~i~ c 
At 


0(Q5/2)^ (20) 


where we have used r = assumed the large 

At limit, and ignored terms of order C>(l/At^). Note 
that Eq. ( |20| is a monotonically increasing function of 
Q, and that this scaling is weak (strictly, 0{Q^) to 
lowest order). Minimizing Eq. ( |20| with respect to C 
gives C = V2(vT^ — 3)/15, and therefore for sufficiently 
small Q-. 


_ 49 + 3v^ 1 _ ^ 

<Pmin ~ (3 yTM)3/2 At At ' 


( 21 ) 


Surprisingly, our small-QST analytics predict that near- 
Heisenberg-limited sensitivities are achievable as Q ap¬ 


proaches zero, with the caveat that r must be very large. 
This is a key qualitative difference between the beam¬ 
splitter and three-mode mixing QST processes. 


B. Numerical Solution 

We complement our small-QST analytic solution - and 
quantitatively investigate the full Q parameter space 
- with numerical phase-space simulations of the three¬ 
mode mixing QST process. Common to all phase-space 
approaches is the conversion of the Heisenberg equations 
of motion to a partial differential equation (PDE) for a 
quasi-probability distribution [3I1|32|. Some approxima¬ 
tion and/or mathematical “trick” is then employed to en¬ 
sure this PDE takes the form of a Fokker-Planck equation 
(FPE) with a positive-definite diffusion matrix, as this 
can be efficiently simulated by a set of stochastic differ¬ 
ential equations (SDEs). Below, we briefly introduce the 
truncated Wigner (TW) and positive-P (P"*") equations 
needed in order to simulate our interferometry scheme 
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(see Sec. IVB1 and Sec. IVB 2[ respectively). The re¬ 
sults of these simulations are presented in Sec. |IV B 3 


1. Truncated Wigner (TW) 

In this paper we have predominantly used the TW 
phase-space method [33HT7] , which is based upon a PDE 
for the Wigner function. Typically, this PDE has third 
and higher-order derivatives that are truncated in order 
to give a FPE; although this is an uncontrolled approx¬ 
imation, it is frequently valid provided the occupation 
per mode is not too small for appreciable time periods 
[551 1^ . Here, the set of TW SDEs corresponding to 
Eqs. Q is 


2. Positive-P (P^) 

In the regime where our low-QST analytics do not ap¬ 
ply, we compared our TW numerical simulations to P'^ 
numerical simulations. The P"*" phase-space method ex¬ 
presses the evolution of the quantum state in terms of the 
Glauber-Sudarshan P representation [3I1I321I1I]. Unlike 
the Wigner representation, the evolution is guaranteed to 
be a FPE provided certain boundary terms are negligible 
(a condition that is valid for sufficiently short times [42 ] ). 
However, a positive-definite diffusion matrix is only en¬ 
sured by doubling the phase-space, which is effected by 
treating the complex amplitudes /? and /3+ = /?*, for ex¬ 
ample, as independent variables. 

A set of P+ SDEs corresponding to Eqs. Q is 


iai=ga 2 P*, (22a) 

id2 = gctiP, (22b) 

iP = ga 2 al, (22c) 

where we have made the correspondences ai(t) —> ai{t) 

and b{t) -T P{t). The initial conditions for these SDEs 
are randomly sampled from the Wigner distribution cor¬ 
responding to the quantum state before the QST process 
(i.e. at time tg) [ID]. Explicitly, Oi is in a coherent state 
of mean number Aft, 0.2 is in a vacuum state, and b is in 
a single-mode squeezed vacuum state. These correspond 
to the initial conditions 


ai{to) = ^/Wt + r]at, (23a) 

0^2{to)=Va2J (23b) 

P{to) = gp coshr -|- rfp sinhr, (23c) 

where rji are complex, independent Gaussian noises sat¬ 
isfying fJ~i=G and r]*r]j = 5^/2. 

The expectation value of any operator function / is 
computed by averaging over solutions to Eqs. (22) with 
initial conditions (23): 


{{f{.a\,al,b\ai,a2,b)}sym) = f {a*i,a*2, P*,ai,a2, P), 

(24) 

where “sym” denotes symmetric ordering [52], and the 
overline denotes the average of simulated trajectories. 
In order to compute the phase sens itivity with the 
information-recycled signal S [Eq. (10)] at the optimal 
operating point, we made use of the following expres¬ 
sions: 


(iV,(U)) = |a,(ti)|2-l/2, (25a) 

(J,(U))=;^, (25b) 

(25c) 

(n(fi)) = -2Im{/3(U)}, (25d) 

= Th^[{p{t^, (25e) 

where i = 1,2 and j = x,z, with J'^ = [ot{{ti)a 2 {ti) + 


«i(^i)« 2 (^i )]/2 and = [\ai{ti)p - \ai\t 2 )\^]/ 2 . 


Si = -igd2P~^ + ^P-igd2Wlit), 

(26a) 

S 2 = —igdiP, 

(26b) 

p = -igd2dp -I- ^J-igwlit), 

(26c) 

dp = igdpP+ ^/igdpW 2 it), 

(26d) 

dp = igdpp'''. 

(26e) 

p+ = igdpdi + y/igw^it). 

(26f) 


where we have made the correspondences apt) -T dpt), 
dl{t) —>■ dp{t), bit) -T Pit), and bpt) -T /3+(t), and 


wpt) are independent complex Wiener noises satisfying 
wpt) = 0 and w*it)wjit') = bijSp — t'). Note that 
Eqs. (26) can be interpreted as either Ito or Stratonovich 
SDEs, since the Stratonovich correction is zero. 

Our initial quantum state corresponds to the initial 
conditions [40] 


aipo) = dpito) = \/f^, (27a) 

02 (^ 0 ) = ap (to) = 0, (27b) 

ppo) = w_(r)ni - i^+ir)n 2 + g, (27c) 

^■^(^o) = -w-(r)ni - v+ir)n2 - g*, (27d) 


where v±ir) = •\/exp(±r) cosh(r)/2, g = (ns -I- 

and Ui are real Gaussian random variables satisfying rU = 

0 and nfrp = Sij. 

Finally, normally-ordered operator expectations corre¬ 
spond to averages over the solutions to Eqs. (26) with 
initial conditions Eqs. (27); for example 


i:fid\,dl, bpai,a2,b):) = fidp, dp, P+,di, d2, P), (28) 

where denotes normal ordering. 


3. Results 


The results of our TW and simulations are plotted 
in Fig. There is excellent agreement between the TW 
and P^simulations in the large-QST regime, and be¬ 
tween the TW simulations and the analytics of Sec. IV A 
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FIG. 4. (Color online) Top panel: Phase sensitivity as a Inac¬ 
tion of QST efficiency for a three-mode mixing QST process 
and an initial condensate of Nt = 10^ atoms. Points joined 
by dotted lines (simply to gnide the eye) are TW simulations, 
whereas the disconnected red diamonds are P^ simulations. 
More specifically, the three sets of square points (teal, orange, 
and green) show TW simulations of the phase sensitivity for 
fixed values of r, and the blue circles show the minimum pos¬ 
sible phase sensitivity [i.e. for optimum r = ropt(Q)] at par¬ 
ticular values of Q. The error bars on the P”*" simulations 
indicate twice the standard error, whereas the standard error 
in the TW simulations is less than the point width. The solid 
magenta curve is the small-QST analytic solution Eq. (20 1 , 
and the horizontal dashed line and dot-dashed line show the 
Heisenberg limit 1 /Nt and minimum possible phase sensitivity 
for a beamsplitter QST process , respectively. 

Bottom panel: Optimal choice of r for a fixed Q as predicted 
by a beamsplitter QST process [r 0 g*(Q) = ln(4W/Q^)/4] 
(black dot-dashed curve), TW simulations (blue circles), and 
small-QST analytics for a three-mode mixing QST process 
^ ln(C^W/Q)l (solid magenta line). 


in the low-QST regime, which give us confidence in the 
validity of our results. 


First, consider the three Ac/) vs Q curves for fixed val¬ 
ues of r. Clearly for r < Tcrit we see that the minimum 
phase sensitivity is always smaller than the minimum Ac/ 
achievable with a beamsplitter QST process. This quan¬ 
titative difference between the beamsplitter and three¬ 
mode mixing QST processes was noted in [23], although 
not explored in any detail. Here, our more comprehen¬ 
sive suite of results show also that for fixed r > rcrit, the 
minimum phase sensitivity occurs at a Q < 1 - a feature 
also observed for the beamsplitter QST results in Fig.|^ 
However, in contrast to the beamsplitter QST process, 
the minimum phase sensitivity decreases as r increases 
(and QST efficiency correspondingly decreases), and in 
fact is close to the absolute minimum « 1.53/iVt once 
the optimal Q < 5%. This is shown very clearly by the 
blue circles in Fig. which are TW simulations of the 
minimum phase sensitivity possible for a hxed Q (i.e. for 
the choice r = ropt(Q) - see bottom panel of Fig. |^. 


C. Optimality with quantum Fisher information 


We now as k th e question: is the information-recycled 
signal S [Eq. (101] the optimal procedure for estimating 
the phase shift (/)? We answer this question using the 
quantum Fisher information if, which places an absolute 
lower bound on the phase sensitivity Acj) > A^qcrb = 
X/^fr^ called the quantum Cramer-Rao bound (QCRB) 
|43l|44]. Crucially, the QCRB applies irrespective of the 
choice of measurement and phase estimation procedure; 
it depends only on the input quantum state. 

In (44] [ 43 , it is shown that when a pure state at 
time ti forms the input to a lossless MZ interferome¬ 
ter, then the quantum Fisher information for estimating 
the relative phase shift between the two arms is simply 
if = AV{Jy{ti)). As discussed in [21], since the initial 
state at time tg is pure, the three-mode mixing QST pro¬ 
cess is unitary, and we are permitting measurements on 
the photons transmitted after the QST process, then the 
state remains pure, and the quantum Fisher information 
for our hybrid atomic-photonic interferometer is the pre¬ 
viously defined J-. Note that this would not be true if 
we only allowed measurements on the two atomic modes; 
then we would have traced over the photon mode, the 
state would be mixed, and this simple expression for the 
quantum Fisher information no longer applies. 

Using Eqs. (15), we obtained an analytic solution for 


the variance in Jy{ti), and therefore the quantum Fisher 
information, valid in the small-r (and therefore the small- 
QST) regime: 


-^C«i(g Q) = Ma r)Nt + A 2 (r, t)N^ + r/N^ + 0 {t^°) 


(29) 





















where 


^4 6 o 

Ai{r, r) = 1 + — [3cosh(4r) + 1] —— [3cosh(2r) + 1] sinh^ r + [35cosh(4r) + 13] sinh^(2r), (30a) 

J. D ^ 0 -L^ 

T-i 

A 2 ir, t) = (e^’’ — l) — ^ 24coshr sinh^ r — 7 cosh(2r) + 3cosh(4r)] 

r® r® 

H-e"" sinh^ r [ sinhr — 10 cosh r + 9sinh(3r) + 6cosh(3r)] H-[8 — 23cosh(2r) + 15cosh(6r)], (30b) 

8 256 

^4 ^8 

Az{r, r) = — (l — sinh^ 32 cosh(2r) + 1] sinh^ r. (30c) 


We also numerically computed the quantum Fisher in¬ 
formation from our TW simulations via Eq. ( |25b ) and 
Eq. (25c) for j = y and J^y = — 

ct2{ti)oil (ti)]/2. 


The phase sensitivity of the information-recycled sig¬ 
nal is compared to the QCRB for three fixed values of 
r in Fig. For Q less than some critical value (ap¬ 
proximately equal to optimal Q), the phase sensitivity 
saturates the QCRB. Close to the optimal Q, where 
the phase sensitivity is minimized, the phase sensitiv¬ 
ity increases, whereas the QCRB is non-increasing with 
increasing Q. However, as shown in Fig. the min¬ 
imum phase sensitivity is never more than a factor of 
~ 1.5 larger than the QCRB. Furthermore, as r increases, 
the minimum QCRB decreases to a minimum of y^/Nt, 
which is only ~ 7% lower than the minimum possible 
sensitivity achievable with the information-recycled sig¬ 
nal Eq. (10). We are therefore justified in calling our 
information-recycling scheme near-optimal] currently an 
optimal phase-estimation scheme eludes us. 


Nevertheless, in the large r regime the optimal phase 
estimation scheme requires a measurement of both the 
atoms at the interferometer output and the transmitted 
photons. Since iV^ = sinh^ r > Nt,hy definition it is im¬ 
possible to achieve complete QST. If only the atoms are 
measured, then the quantum Fisher information depends 
on the reduced density matrix of the atoms alone. Here 
the quantum Fisher information is not simple to compute 
(it requires the diagonalization of the symmetric logarith¬ 
mic derivative), although it is guaranteed to be less than 
or equal to AV{Jy{ti)) [ 33 ]. Furthermore, since imperfect 
QST acts, to a first approximation, as a linear loss mech¬ 
anism, Heisenberg scaling rapidly reverts to 1/^/Nl for 
even slight departures from Q = 1 [H] . This is consistent 
with the results of [23] . which showed that without in¬ 
formation recycling imperfect QST severely reduces the 
enhancement due to the squeezed light. Consequently, al¬ 
though information recycling is not strictly optimal, the 
optimal phase estimation procedure must incorporate in¬ 
formation from both atomic and photonic measurements 
- as information recycling does - in regimes where only 
Q < 1 is possible. 



FIG. 5. (Color online) Comparison of the QCRB and 
the phase sensitivity of the information-recycled signal for a 
three-mode mixing QST process and an initial condensate of 
Nt = 10^ atoms. Points indicate the results o f TW simula¬ 
tions; squares are A(f> calculated using S [Eq. ([l0|], and tri¬ 
angles (joined by a dashed line to guide the eye) indicate the 
QCRB. The solid magent a an d blue curves indicate t he s mall- 
QST anal 3 rtic curves Eq. (191 and l/y/iFa^i [see Eq. |2^], re¬ 
spectively, for r = 6.31. Here, there is a very good agreement 
between the TW simulations for the QCRB and the small- 
QST analytic solution, although there is a slight deviation 
between the TW simulation for the phase sensitivity of the 
information-recycled signal and Eq. ( |19| ). Since the standard 
error of all points is less than the point width, this discrep¬ 
ancy is caused by slight deviations from a physical Wigner 
functiorQ due to the uncontrolled approximation associated 
with the TW simulation method. 


^ A breakdown in the pure state expression for the quantum 
Fisher information, tF = 4V{jy{ti)), is not possible; since the 
QST dynamics are unitary, the Wigner function is guaranteed 
to remain pure, even under the TW approximation [46]. 


V. CONCLUSION 


In most quantum technological applications that in¬ 
corporate QST between two quantum systems, a larger 
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FIG. 6. (Color online) Truncated Wigner simulations of the 
minimum phase sensitivity A(j>^in{r) = miriQ A4){r, Q) and 
minimum QCRB for a three-mode mixing QST process and 
initial condensate of Nt = 10^ atoms. As r gets large, the 
QCRB asymptotes to \/2jNt, which is only ~ 7% lower than 
the minimum possible sensitivity « 1.53/At achiev¬ 

able with the information-recycled signal [see Eq. ( |21[ )]. The 
standard error in each simulation is less than the point width; 
consequently, the slight “wiggles” in the magenta points at 
large r are systematic errors due to a breakdown in the uncon¬ 
trolled approximation underlying the TW simulation method. 


xlQ- 


* QCRB 

• Information-recycled signal 



QST efficiency is seen as more desirable. In this paper, we 
have provided a clear counter-example to this intuition. 
Specifically, we showed that if a hybrid atomic-photonic 
interferometer is enhanced by single-mode squeezed op¬ 
tical vacuum and information recycling (see Fig. [^, and 
the information-recycled signal Eq. ( (l0| is constructed, 
then for moderate to large levels of squeezing the best 


phase sensitivity requires an imperfect QST efficiency. 

If the QST process is a beamsplitter, then this mini¬ 
mum phase sensitivity is , for total atom 

number Nt. As shown in Fig.[^ this minimum is obtained 
provided r = > Ccrit., where rgg* « ln(4At/Q^)/4 

and Tcrit = ln(4iV()/4. In contrast, we showed that a 
three-mode mixing QST process yields a phase sensitiv¬ 
ity close to the Heisenberg limit; counterintuitively, this 
is only achievable in the very small-QST, large squeez¬ 
ing regime (see Fig. |^. We compared this result to the 
QCRB, and showed that our scheme is nearly optimal 
(see Fig. and Fig. [^. In particular, the QCRB pre¬ 
dicts a minimum phase sensitivity of V^/Nt, which is 
only ~ 7% lower than the minimum achievable with our 
information-recycled signal. 

Unfortunately, the squeezing parameters required in 
order to realize these three-mode mixing QST dynam¬ 
ics are much larger than those achievable with current 
technology. We therefore see little prospect of an exper¬ 
imental observation of this effect in the near term. Nev¬ 
ertheless, the results of this paper are still of theoretical 
interest, as they demonstrate that information recycling 
is capable of exploiting metrologically useful correlations 
in regimes where quantum metrology would naively be 
deemed impossible. 
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